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Abstract 

Filter bank-based multicarrier communications (FBMC) have recently 
attracted increased interest in both wired (e.g., xDSL, PLC) and wireless 
(e.g., cognitive radio) applications, due to their enhanced flexibility, higher 
spectral efficiency, and better spectral containment compared to conven- 
tional OFDM. A particular type of FBMC, the so-called FBMC/OQAM or 
OFDM/OQAM system, consisting of pulse shaped OFDM carrying offset 
QAM (OQAM) symbols, has received increasing attention due to, among 
other features, its higher spectral efficiency and implementation simplic- 
ity. It suffers, however, from an imaginary inter-carrier/inter-symbol in- 
terference that complicates signal processing tasks such as channel esti- 
mation. This paper focuses on channel estimation for OFDM/OQAM 
systems based on a known preamble. A review of the existing pream- 
ble structures and associated channel estimation methods is given, for 
both single- (SISO) and multiple-antenna (MIMO) systems. The vari- 
ous preambles are compared via simulations in both mildly and highly 
frequency selective channels. 
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1 Introduction 



Orthogonal frequency division multiplexing (OFDM) has become quite popular 
in both wired and wireless communications [71 140[ 186] . mainly because of its 
immunity to the frequency selectivity of the channel, which allows a significant 
increase in the transmission rate [531 EH]- Using the cyclic prefix (CP) as a 
guard interval against intersymbol interference (ISI) , OFDM manages to turn a 
frequency selective channel into a set of parallel flat channels with independent 
noises. This greatly simplifles both the estimation of the channel and the re- 
covery of the transmitted data at the receiver. However, these advantages come 
at the cost of an increased sensitivity to frequency offset and Doppler spread. 
This is due to the fact that, although the subcarrier functions are perfectly time 
limited, they suffer from spectral leakage in the frequency domain and hence 
inter-carrier interference (ICI) results under nonideal conditions. Moreover, the 
inclusion of the CP entails a waste in transmitted power as well as in spectral 
efficiency, which, in practical systems, can go up to 25% '7^. 

An alternative to CP-OFDM, that can mitigate these drawbacks, is given by 
filter bank-based multicarrier (FBMC) systems, which have recently attracted 
increased interest [73] for both wired (e.g., power hue comms (PLC) [TOllTSl fT]) 
and wireless (e.g., cognitive radio [SnilHSlIill and DVB-T [3]) applications. For 
an excellent presentation of FBMC systems, including a review of recent appli- 
cation areas and related standardization activities, the reader could refer to [31) . 
FBMC enjoys enhanced flexibility, higher spectral efficiency, and better spec- 
tral containment compared to CP-OFDM [31(51]. Special attention has been 
recently paid to an OFDM filter bank-based variant employing offset quadrature 
amplitude modulation (OQAM), known as OFDM/OQAM [5511211111117]. This 
was, in its principle, first proposed back in the 60's [201175], but the interest in 
this scheme was only recently revived due to its desirable features including the 
potential for maximum spectral efficiency as well as implementation simplicity 
over alternative FBMC schemes [121 [T51 [75] . It consists of pulse shaping, im- 
plemented via an IFFT/FFT-based efficient filter bank, and the subcarriers are 
loaded with staggered OQAM symbols, i.e., real symbols at twice the symbol 
rate of OFDM/QAM [SgHHO]. This allows the pulses to be well localized in both 
the time and the frequency domains |18) . thus increasing the system's robust- 
ness to frequency offsets [21 177] and Doppler effects [IHl jZH [571 |3B] and provid- 
ing better spectral containment in bandwidth sensitive applications [51 [HI HI] ■ 
Furthermore, although the OQAM-based OFDM scheme can be seen to ex- 
hibit similar peak-to-average power ratio (PAPR) performance with classical 
OFDM/QAM [55], the presence of spectral leakage in the latter may ultimately 
generate higher peak power values [5S]. Moreover, because of the good time 
localization, a CP does not have to be included in the OFDM/OQAM trans- 
missionQ which leads to still higher transmission rates [501 El [71] . 

However, the subcarrier functions are now only orthogonal in the real field, 
which means that there is always an intrinsic imaginary interference among 
(neighboring) subcarriers and symbols [45) . This makes the channel estima- 
tion task for OFDM/OQAM systems more challenging. A number of training 
schemes and associated estimation methods have recently appeared in the lit- 
erature, including both preamble-based [151[Sl[nSl[lIl[ll[lSl[Sl[S3]and scat- 

^ Nevertheless, this advantage was partly given up in I66| and a CP-based OFDM/OQAM 
system was proposed for the sake of facilitating the data reception process. See also 1211 . 
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tered pilots-based gSllIlllHllMllllSllHilMllHllSIllSTlISS] ones. These can 
be roughly characterized as aiming at cancelling/avoiding the undesired inter- 
ference or constructively exploiting it to improve the estimation performance. 
The focus of this paper is on OFDM/OQAM channel estimation based on a 
known preamble. Existing preamble structures and their associated estimation 
methods are reviewed, for both the single-antenna (SISO) and multiple-antenna 
(MIMO) scenarios. The problem of preamble optimization to achieve minimum 
estimation error is also shortly discussed, based on recently reported results. 
The estimation performances of the principal methods are compared via simu- 
lations for both mildly and highly frequency selective channels. 

The rest of the paper is organized as follows: In Section [2l the discrete- 
time baseband equivalent model for the OFDM/OQAM system is described. 
Section |3] reviews the main approaches proposed in the literature for structuring 
the preamble and exploiting it for SISO channel estimation. The MIMO case is 
addressed in Sectional Comparative simulation results are reported in SectionjS] 
Section [6] concludes the paper. 

Notation. The complex conjugate of a complex number z is denoted by z*. 
Also, J = v^— 1. 3?(z) or z^ denotes the real part of z, while its imaginary part 
is written as ^(z) or z^. The set of real numbers is denoted by IR. Boldface 
lower and upper case letters are used for vectors and matrices, respectively. The 
nth-order identity matrix is denoted by /„. || • || is the Frobenius norm. The 
expectation is denoted by E{-). is the (left) Kroncckcr product. 

2 System Model 

The (QAM or OQAM based) OFDM modulator output is transmitted through 
a channel of length Lh, which is, as usual in block transmissions, assumed to 
be invariant in the duration of an OFDM symbol. At the receiver front-end, 
noise is added, which is assumed Gaussian with zero mean and variance a^. For 
CP-OFDM, the classical configuration is assumed |H3]- The discrete-time signal 
at the output of an OFDM/OQAM synthesis filter bank (SFB) is given by [SDpl 

A/-1 

S(0 = X! ^^^n,ng7n,n{l), (1) 
m—O n 

where dm,n are real OQAM symbols, and 

with g being the real symmetric prototype filter impulse response (assumed 
here of unit energy) of length Lg, M being the even number of subcarriers, and 
'fim,n = <^o + f + n) mod tt, where ipa can be arbitrarily chosen |80| . In this 
paper, ipm,n is defined as {m + n)^ — mmr as in [SD]. The filter g is usually 
designed to have length Lg = KM, with K being the overlapping factor. The 
double subscript (•)m,n denotes the (m,n)-th frequency-time (FT) point. Thus, 
m is the subcarrier index and n the OFDM/OQAM symbol time index. 

^Alternative configurations exist, generating a real OFDM/OQAM signal; see I39II84|[67] . 
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The pulse g is designed so that the associated subcarrier functions g,n,7i are 
orthogonal in the real field [TH] , that is 



(Offp,g(0 { = S„i^pSn,q, (2) 



where Sij is the Kronecker delta. This implies that even in the absence of chan- 
nel distortion and noise, and with perfect time and frequency synchronization, 
there will be some intercarrier (and/or intersymbol) interference at the output 
of the analysis filter bank (AFB), which is purely imaginary (the notation in 
[62] is adopted): 

T.9mAl)9;,gil)^j{9rra% (3) 
I 

and known as intrinsic interference [45) . Making the common assumption that 
the channel is (approximately) frequency flat at each subcarrier and constant 
over the duration of the prototype filter [S2] , which is trueH for practical values 
of Lh and Lg and for well localized 5's, one can express the AFB output at the 
pth subcarrier and qth OFDM/OQAM symbol as: 



where flp,g is the channel frequency response (CFR) at that FT point, and Ip,q 
and rjp^q are the associated interference and noise components, respectively, rjp ^ 
has been shown to be stationary (in fact, also Gaussian with zero mean and 
variance a^) and correlated among subcarriers (see [H]). In the sequel, this 
correlation will be mostly neglected, as it is this case with well localized AFBs. 

For pulses 9 that are well localized in both time and frequency, the interfer- 
ence from FT points outside a neighborhood flp^q of (p, q) (excluding {p, q) itself) 
is negligible. If, moreover, the CFR is almost constant over this neighborhood, 
one can write ^ as 

Up^q ~ Hp qCp q -\- rjp q (5) 



where 



Cp.g — dp q -|- ] ^ ^ dm.n{9)m'^,n ~ dp,q + JUp,q (6) 
(m,n)gf2p^. 



is the virtual transmitted symbol at (p, q) , with 



(m,n)Gf!p ,j 



^This is not an accurate subchannel model in case M, the number of subcarriers, is not 
large enough with respect to the frequency selectivity of the channel |70| (see Section [5] for 
such examples). This is also the case in high mobility scenarios, where M should be small, to 
cope with ICI from frequency dispersion 1321 1331 l34l . 
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being the imaginary part of the interference from the neigboring FT points. 
When known pilots are transmitted at that FT point and its neighborhood 
Op^g, the quantity in ^ can be approximated. This can then serve as a pseudo- 
pilot to compute an estimate of the CFR at the corresponding FT point, 
as, for example, 

yp,q ^ jj I ''lp,q 



Hp^q — 



(8) 

This observation underlies the so-called Interference Approximation Method 
(lAM) to be desribed in detail in the next section. 

The most common assumption is that, with a well time-frequency localized 
pulse, contributions to Ip^q only come from the first-order neighborhood of {p, q), 
namely ilp^q = {{p ±l,q± 1), {p, q ± 1), {p ± 1, q)} (see Fig. [1]). A special case 
is given by ^p^q — flp ^ — {{p ± 1, q)}. This arises when we place three adjacent 
pilot tones at positions (p — I, q), (p, q), (p + l,q) and zeros at the rest of the 
first-order neighborhood positions or when we place nonzero pilot tones at all 
positions in the training vector and zero vectors around it. 

If the neighbors of (p, q) carry unknown (data) symbols, one cannot (directly) 
approximate the imaginary interference in ([B]). However, by properly choosing 
one of the neighboring symbols, say at the point (r, s), this interference can be 
forced to zero. Then the pseudo-pilot in ([6]) becomes real, equal to dp^q. The 
pilot at (r, s) is then known as a help pilot [3S] 

The above formulation can be easily extended to the MIMO case. Consider 
an A't X iVr MIMO system, with identical SFBs ( AFBs) at each transmit (receive) 
antenna (see Fig. [5]) . Then, adopting the same assumptions as before, one can 
write an equation analogous to ([5]) for each receive antenna j = 1,2, . . . , N,-: 



(9) 



where 77^^* is the Af-point CFR of the channel from the ith transmit antenna 
to the jth receive antenna, ^ is the corresponding virtual symbol, namely 



dl.q+J E 

(m,n)eOp ,j 



m.n \ci / m,n 



J^p,q' 



and 7]^ q denotes the corresponding noise component. Clearly, the latter is un- 
correlated among different receive antennas. Collecting the above equations for 
all receive antennas, results, for each FT point, in an input-output equation 
similar to that for MIMO-OFDM [521: 



where 





yp,q ~ H 


p,q'^p,q 


^p,ql 


/ 1 

p,q ^ 


^pk 

p,q 


rrl.2 

^P.q 

TT2.2 

p.q 


■ ■ ■ H^'^^ 

■ ■ ■ 

p,q 




rrJV,.l 
L ^P,q 


^P,q 


. . . frN,,Nt 

p.q 



(10) 



^This idea was later proposed again in 1481 . 
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is the MIMO CFR at that point and 
with 

d - r rfi (P • • • d^' 1 ^ 

V^l ~ L P,9 P,9 P,<3 J 

and itp ^ defined similarly. 



3 Preamble-based Channel Estimation Methods: 
The Single- Antenna Case 

Known (training) symbols for the purpose of channel estimation can be gath- 
ered in the beginning of a frame and/or placed at isolated FT points throughout 
the frame. These are the two principal pilot configurations, known as preamble 
and scattered pilots, respectively. The focus of this paper is on the design of 
the preamble!! for OFDM/OQAM systems and its use in estimating the chan- 
nel. The difficulty in this task, compared to OFDM/QAM, lies in estimating 
a complex CFR in the presence of the imaginary intrinsic interference among 
neighboring FT points. The preamble design and the associated estimation 
method should take this scenario into account in some appropriate way. Three 
different approaches for preamble-based channel estimation, including their vari- 
ants, are reviewed in this section, for SISO systems. The first approach. Pairs of 
Pilots (POP), relies on simple algebraic relations for the input/output samples 
in a number of (in practice consecutive) time instants. The second one, In- 
terference Approximation Method (lAM), aims at approximating the intrinsic 
imaginary interference from neighboring pilots and hence constructing complex 
pseudo-pilots as above, to accommodate the complex CFR. A number of vari- 
ants of I AM are considered. In this context, recent results on designing optimal 
- in the sense of minimizing the channel estimation error - preambles are also 
recalled. An alternative approach, that of cancelling or avoiding interference, is 
then discussed. The extension of some of these methods to the MIMO case is 
addressed in the following section. 

Note that the methods presented here can be applied to any OFDM/OQAM 
filter bank, provided, of course, that the prototype filter is sufficiently well local- 
ized in time and frequency and the number of subcarriers is large enough that 
the models in eqs. ([5]), ([TU|) hold with sufficient accuracy. Thus, methods aim- 
ing at enhancing the channel estimation performance through an interference- 
minimizing filter design are omitted; see, e.g., |22j[68j. Moreover, we do not 
include methods based on preambles that are not completely known to the re- 
ceiver (as in, e.g., |35j. where the training sequence is superimposed on the data) 
or are too spectrally inefficient (e.g., made of isolated pilots surrounded by a 
large neighborhood of zeros, as in [7T1[711[73]) or especially suited to wideband 
sub-channels, being significantly more complicated (e.g., [3 3) ) or significantly 
longer than those considered here (e.g., ^38 ^ 137 ] ). 

In the sequel, and in addition to the assumptions stated in the previous 
section, the channel is assumed to be time invariant during the preamble period. 

^Same principles apply if the pilot blocks arc placed in the middle of the frame (midambles). 
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3.1 Pairs of Pilots (POP) Method 

This method, proposed in [52], aims at computing a channel estimate by us- 
ing ^ at two different (in practice consecutive) time instants, 91,(72, to con- 
struct a system of equations for the real and imaginary parts of Hp q. Here, it 
is described in an alternative, equivalent way |23| . 

With the noise being neglected, denote by Wp^q = 1/Hp_q the corresponding 
zero forcing (ZF) equalizer coefficients. Then, one can write 



yp,qi ^p,q 



- JUp,c 



and therefore 



yp,qi " p,qi 
yp,q2 P,<32 



yp,qi " p,qi 
yp,q2 P,q2 



= d 



'p,qi 



'P,q2 



Noting that Wp,q-^ ~ W^p,g2 '^^^ to the (almost) time invariance of the CFR, 



yp,qi 



y^,q 



p,qi 
p-m 



J'P,<32 



hence 



p,qi 



Vp,qiyp,q2 yp,qiyp,q2 



This can be more compactly written as 



yp,qi'^P,q2 ~ 2/p,q2"'P.9l 



yp,qi'^P-q2 yp,q2"'P'<ll 
rj — j;R rj 



p,qi 



qi 



^p-qi Vp.t 



^P,q2yp,qi 



^(yp,gi2/p,i32) 



(11) 



In practice, the preamble would comprise the first two OFDM/OQAM symbols, 
that is, qi — and q2 — I- In that case, and with the preamble suggested 
in [52], where dp^ ~ (—1)^ and the second symbol is all zeros, the above would 
simplify to 

i-Wp,i 



P,Q =J 



^(2/p,o2/p,i) 



(12) 



The CFR would then be computed as Hp^o = l/Wp^. 

An advantage of the POP scheme, besides its simplicity, is that it does not 
explicitly depend on the employed prototype filter. However, it must be em- 
phasized that the above derivation only holds when noise is negligible. One can 
see that (see [52]), in the presence of noise, the method can have unpredictable 
performance since the degree of the noise enhancement in general also depends 
on unknown (hence uncontrollable) data. 



3.2 Interference Approximation Method (lAM) 

If the interference jUp,q in ([5]) is only due to the immediate neighbors of (p, q) and 
these FT points carry training (hence known) symbols, then one can compute 
an approximation of this interference, as pointed out above. This can then be 
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used to construct the complex pseudo-pilots Cp^q and subsequently use them to 
get a channel estimate in exactly the same way as in CP-OFDM (see ([5])). This 
method is called (for obvious reasons) the Interference Approximation Method 
(lAM) and requires that all input symbols in the immediate neighborhood of 
{p,q) be known. This implies that, for each subcarrier, (at least) 1.5 complex 
symbols are required for training, representing an increased (though not signif- 
icantly) overhead compared to CP-OFDM. The training would need to extend 
over more than three time instants for less well time-localized filters. 
Writing (O out. 



Vp, 



9 



one can see that the pilots should be so chosen as to minimize the effect of 
the noise term. In [611 [62] , it is proposed that these be chosen so as to have 
pseudo-pilots of maximum magnitude. To this end, the surrounding training 
symbols should be such that all terms in ([7]) have the same sign so that they 
add together. Moreover, this should happen at all frequencies p. 

To this end, one needs to know the interference weights for the neigh- 

bors {m,n) G ftp^q of each FT point (p, q) of interest. These can be a priori 
computed based on the employed prototype filter g. Moreover, it can be shown 
(see the Appendix) that, for any choice of these weights follow a specific pat- 
tern, which, for the definition of ipo adopted here, and for all q, can be written 
afH 

{-l)Pe -(-l)^e 

H-lYl dp^q (-1)^7 (13) 

{-Ife -(-l)^e 

with the horizontal direction corresponding to time and the vertical one to 
frequency, as in Fig. [TJ The above quantities can be shown (see the Appendix) 
to be given by 



; 277 



L„-l 



-J M 2 

1=0 

i„-l 



E 5'(0e^'^' (14) 



7 = 



E -9(0.9 ('-V) (15) 



i_M_ 



2 , 



6 = -je-'T^^ E .9(0.9 ^) e'^' (16) 



l=M/2 
L„-l 



J2 9(1)9 (l~^)e^^'^' (17) 

l=M/2 



°This is without loss of generality as analogous patterns result for other definitions of ipo 
(see, e.g., [26]). 
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and are (absolutely) smaller than one, with /3, 7, 5 > 0. Generally, /?, 7 > (5 > |e|. 
For example, 7 = 0.553, (3 = 0.25, S = 0.2172, and e « 0.0004 for the g's 
designed in 11 with M = 512, = 3. 

3.2.1 The lAM-R Method 

In order to simplify the task of generating pseudo-pilots of large magnitude, it 
was suggested in [SlllSllSn] to place nulls at the first and third OFDM/OQAM 
symbols of the preamble, that is, dp^ — dp^2 = 0, p = 0, 1, . . . , M — 1. Then, 
the imaginary part of the pseudo-pilot at (p, 1) will only come from the symbols 
at the positions (p ± 1,1). Obviously, these pilots must be OQAM symbols of 
maximum modulus, d. Moreover, in view of the pattern (I13p . they should sat- 
isfy the condition (ip+1,1 = — rfp-i,! for all p, in order to yield pseudo-pilots of 
maximum magnitude, namely \dp,i + j2l3dp^i^i\ = dy'l + AfP. This preamble 
structure is shown in Fig. ^a.), for the example of M = 8 and with OQPSK 
modulation. This channel estimation method is called IAM2 in |S2], to distin- 
guish it from lAMl, where all symbols, at all three time instants, are randomly 
chosen. To prevent possible increase of the peak to average power ratio (PAPR) 
with IAM2, a random choice of the middle symbols in IAM2 is also proposed 
in [55], giving rise to the IAM3 method. In the sequel, we will focus on IAM2, 
due to its ability to maximize the magnitude of the pseudo-pilots with the aid 
of real pilot symbols, and use the name lAM-R, coined in [65] to signify the fact 
that its pilot symbols are real. 

3.2.2 The lAM-I Method 

It can be verified from ^ that one can do better in maximizing the magnitude 
of the pseudo-pilots if imaginary pilot symbols jdp.i of maximum modulus d are 
also allowed in the preamble. The reason is that the corresponding pseudo-pilot 
would then be imaginary and hence of larger magnitude than in lAM-R, namely 
\dp,i + Up^i\, provided that the signs of the pilot symbol and its neighbors are 
properly chosen so that > 0. Such an lAM scheme was first proposed 

in [5S], and was called lAM-I to signify the presence of imaginary pilots. The 
middle preamble vector consists of triplets, with each one of them following the 
above principle, and the pilots in each triplet are otherwise selected in random 
and independently of the other triplets. Hence, imaginary pseudo-pilots result 
in only one third of the subcarriers, with magnitude d{l + 2/?), whereas the rest 
of them deliver complex pseudo-pilots of smaller magnitude, d\{l + /?) + An 
example, for M = 8 and OQPSK symbols, is depicted in Fig. [2I^b). Note that 
the above preamble is, strictly speaking, not OQAM. Nevertheless, it can still 
be fed to the synthesis filter bank and perfectly reconstructed at the analysis 
bank. 

3.2.3 The lAM-C Method 

With the preamble proposed in |65j , imaginary pseudo-pilots will result in only 
one third of the subcarriers used, whereas the rest of the subcarriers will deliver 
complex pseudo-pilots. With a slight modification, one can improve upon this 
and obtain pseudo-pilots that are either purely real or imaginary at all the 
subcarriers. The idea is to simply set the middle OFDM/OQAM symbol equal 
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to that in lAM-R but with the pilots at the odd subcarriers multiphed by j, 
as shown in [53] and independently in [H]. Again, for the M = 8 OQPSK 
example, this will result in the preamble shown in Fig. |3t^c). This preamble was 
called lAM-C in 

3.2.4 Optimal Preambles 

In [55], the lAM-C preamble was derived as a solution to the problem of optimal 
3-symbol preamble design with zero guard symbols, where optimal is in the sense 
of generating pseudo-pilots of maximum magnitude. Optimality of preambles 
of this kind, in the sense of attaining the minimum possible mean squared error 
(MSB) with a given transmit energy budget, was recently studied in |49j and 
the results were later complemented in [50'. Two different configurations were 
considered: sparse, where only isolated subcarriers carry pilots, with the rest 
of them being nulled, and full, with all frequencies being occupied with nonzero 
pilots. It was proved that: 

1. The pilot-carrying subcarriers in an optimal sparse preamble are equidis- 
tant with the pilots being equipowered. 

2. With (prn,n defined as above and including the phase factors e-''^™-^ in the 
pilot symbols, they should be all equal but with alternating signs in an 
optimal full preamble^ 

One can readily see that the lAM-C preamble - proved to be optimal in [55] - 
follows the general structure of an optimal full preamble of this kind, as it was 
derived in [49 , ,50 , . For example, multiplied with e-''^™-^ the middle vector in 
Fig. |3i;c) becomes j, -j, 

3.2.5 The Extended lAM-C (E-IAM-C) Method 

The optimality of lAM-C is restricted to the 3-symbol preambles with nulls at 
their sides. One could wonder whether an alternative preamble structure which 
employs the side symbols as well could yield pseudo-pilots that are stronger 
than those of lAM-C. It turns out that such a preamble can indeed be formed 
in an analogous way as previously, by employing the left- and right-hand sides 
of the neighborhood as well [52^1^ Thus, at an odd-indexed subcarrier p, 
with the pilot ±dj in the middle, one should place at the right hand side, 
and its negative, ±d at the left hand side. These side pilots will then contribute 
±d7j each to the interference for the central pilot. In an analogous way, at an 
even-indexed subcarrier p, the middle pilot, say ±d, will get a total interference 
of ±2djj from its {p, 0) and {p, 2) neighbors, if these are chosen as ±dj and 
=F(ij, respectively. It can then be verified that, because of the interference 
components with weights ±5 cancel each other, and hence have no contribution 
to the pseudo-pilot. However, this loss is not significant since S is significantly 
smaller than 7. The resulting pseudo-pilots are again either purely real or 

^With ipm.n = + n)-K jl instead, these symbols must be all equal. See 511 for more on 
this, independent of the ipm,n definition. 

*This idea seems to have first been investigated in 1261 Section 5.4]. It must be noted 
however, that, for a number of pulse types (including those of 1111 employed for the simulations 
in this paper, where 7 > /3 > <5 and 5 < 0.5), the preamble of I52| results in pseudo-pilots of 
a larger magnitude than in the preambles developed in | 26| . 
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imaginary, with magnitude c?| 1 + 2(/3 + 7 + 2e) | , which is clearly larger than that 
achieved by lAM-C. For the example filters employed here [TT], with M = 512 
and K ~ 3, the corresponding magnitudes are 2.6076d and 1.5d. An example 
of the preamble configuration for M = 8 and OQPSK modulation is given in 
Fig. ^d). One can see that the left hand column is built by repeating each 
quadruple of the middle pilots in reverse order, while the right hand column is 
simply the negative of the left hand one. 

The above hold for prototype filters having e > 0, as it is most often the 
case. Whenever e < 0, as it is the case for Extended Gaussian Function (EGF)- 
based filter banks for example [60], the E-IAM-C preamle can be modified to 
yield pseudo-pilots of magnitude •\/(l + 2/3)2 _|_ 4^^ „ 2e)2. Fig.jSje) shows such 
a preamble configuration for M — 8, with OQPSK modulation. 

It must be emphasized here that, since all three OFDM/OQAM symbols of 
the E-IAM-C preamble are nonzero, it will require more power to be transmitted 
than the previous ones. This fact is fairly taken into account in the simulation 
experiments (see Section [5]), along with possible interference coming from the 
data part of the frame. 

3.3 Interference Cancellation/ Avoidance Methods 

An alternative, more straightforward approach to cope with the intrinsic inter- 
ference is to so design the preamble and the associated estimation method as to 
somehow cancel or avoid it. This is in contrast to the lAM approach presented 
above, where interference is taken advantage of to improve the channel estima- 
tion accuracy. The methods following this alternative approach can be roughly 
categorized as follows. 

A direct way to cancel the interference is to simply null the data surround- 
ing the pilot FT point. Note that this idea was already used in some of the 
lAM preambles above, where guard OFDM/OQAM symbols were inserted to 
(approximately) null the interference from preceding and following data. The 
problem can be further simplified if zeros are also transmitted at the middle 
symbol, at all odd- (or all even-) indexed subcarriers. In this way, the CFR can 
be estimated at the even-(odd-)indexed subcarriers just like in CP-OFDM. The 
channel at the rest of the frequencies will then have to be estimated via inter- 
polation |36j . Such a preamble design was proposed in |42) and an example (for 
M — 8) is given in Fig.SJa). The same idea was essentially used in [751 [73], and 
it was also independently included in the preamble evaluation study of [24j (for 
MIMO). In the latter, only every third subcarrier was modulated, in accordance 
with the pilot configuration suggested in WiMAX DL-PUSC with segmentation 
[7]. A less direct design, which relies on the symmetries in (fT3)) to cancel the 
interference from adjacent subcarriers, was proposed in |48[ and is depicted 
in Fig. UKb) for the M — 8 example. 

A more spectrally efficient preamble design method, which also takes ad- 
vantage of the symmetries in the neighboring times in p3l) to transmit data 
at these FT positions instead of nulls, was suggested in [5S]. The way this is 
achieved is by considering only the first-order neighbors and canceling the in- 
terferences coming from them in pairs. Thus, following the pattern in (|13p. the 
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data modulating the first-order neighbors would be structured as follows: 

a b c 

d X d 
—a b — c 

Note that the idea is similar to that of the help pilot |45J, although, here, only 
four, instead of seven information symbols can be transmitted per pilot. In fact, 
this idea is proposed in |88| in the context of scattered pilots. One can however 
easily use it to design a preamble with the property of interference cancellation, 
by simply applying the above pattern for every third subcarrier and resorting 
to interpolation for the rest of them. This would result in a preamble contain- 
ing about data symbols; this is a considerable gain in spectral efficiency 
compared to the methods previously described. To obtain better channel esti- 
mates at the intermediate frequencies, one could even apply this in every second 
subcarrier. This is made possible by the fact that, as shown in (|13p. the interfer- 
ence weights are the same for every second subcarrier. However, such a preamble 
construction would impose some additional constraints to the data that can be 
transmitted, reducing the data symbols to only 3 -I- -y. Fig. |4l^c) provides an 
example for M = 8 and pilots placed at the even-indexed subcarriers. 

Preambles of an even higher spectral efficiency can be conceived, comprising 
only one reference OFDM/OQAM symbol surrounded by completely unknown 
data. Hence, no guard symbols or data of a specific structure are transmitted 
along with the known symbols as in the previous cases. This, however, implies 
that the interference from surrounding symbols is completely unknown and its 
elimination becomes an important issue. The approach proposed in |41j is to 
follow an iterative procedure, where in each iteration the available channel esti- 
mate is used to approximate the surrounding data and vice versa. The original 
channel estimate is only based on the preamble symbol, ignoring the interfer- 
ence effect. The following channel estimates are computed on the basis of both 
pilot symbols and estimated data symbols (e.g., using least squares (LS)). The 
channel and data estimates are improved from one iteration to the other. In the 
simulation results reported in [41^, the iterative approach is shown to achieve 
better performance than the method relying on guard symbols, in only two it- 
erations. More iterations are shown to lead to further improvements. Such an 
iterative approach was also proposed in albeit for isolated (scattered) pi- 
lots. Observe that one can so choose the pilot symbols so as to, at least, cancel 
the interference from neighboring pilots, by respecting the symmetry of (|13|) . 
An example, for M = 8, is given in Fig. HJd). A similar preamble is used in 
the method of [25], where the channel and data symbol estimates are used to 
approximate the interference and subtract it from the AFB output. 

4 Preamble-based Channel Estimation Methods: 
The Multiple- Antenna Case 

4.1 MIMO POP Method 

The POP method of Section 13.11 can be easily extended to the MIMO case, 
as follows. Consider (|10l) with noise being neglected and assume ^ is left- 
invertible for all p,q (hence N^- > Nt). Denote by W the (per-subcarrier) 
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ZF equalization matrix, that is, = -ffp^g, where H^^^^ is the left inverse 

of Hp g. Write ([TU| in the form W^^ypg — c^^^ for 2N^ time instants, say 
0' = 0, 1, . . . , 2N,: — 1, assuming that the channel (and hence W) can be viewed 
as invariant in this time period. Then 

11^ — 11^ — d 

pfiUp.q p.OHp.q ~ "'p..q^ 

for g = 0, 1, . . . , 2N^ — 1. Collecting these equations into one, one can write 



-yl.Q 



-vis 



yp,2N,-i 

^yp,2N,-l 



■p,0 



d 



'p,2N^-l 



or 



The ZF equalizer can then be computed as 



W- 



DpY- 



p,0 '^^ p,0 

and hence the channel matrix can be found as the right inverse of 

1 



J 



Remark. It must be noted that this scheme is no longer associated to the concept 
of "pairs of pilots" and the name "POP" can only be used here as a reminder 
of its SISO counterpart. 



4.2 MIMO lAM 

As in MIMO-OFDM, one can consider three different patterns for preamble 
training (see, e.g., [71 Section 5.7]): independent (one antenna at a time), scat- 
tered (transmitting at different frequencies), and orthogonal (transmitting or- 
thogonal signals at the same times and frequencies). lAM preambles of the 
latter category can be easily constructed from their SISO counterparts if a pat- 
tern analogous to that used for MIMO-OFDM in ^82j is followed. Such an 
extension of lAM to the MIMO setup was reported and analyzed in [53], and is, 
in fact, the first attempt to address full preamble-based channel estimation in 
the MIMO-OFDM/OQAM context. A sparse preamble design was first consid- 
ered in [21] and was based on earlier work on orthogonal MSE-optimal training 
design for MIMO-OFDM systems. This will be described in the next subsection. 

It is clear from dTO]) that one wiU need at least Nt (nonzero) OFDM/OQAM 
symbols to estimate the CFR matrix. Assume the channel is time invariant 
during a period of 2Nt + 1 OFDM/OQAM symbols (SiVt for E-IAM-C). Take 
the example of A^t = 2, for the sake of simplicity. One antenna can use the same 
preamble as the one used in the SISO case, but with a repetition, as shown in 
the example of Fig.jSfa), for lAM-C (similarly for the other lAM variants). The 
other antenna uses the same preamble, but with reversed signs at the second 
nonzero OFDM/OQAM symbol (cf. Fig. [Sib)). Then, writing ^ at times 
g = 1, 3 (g = 1, 4 for E-IAM-C) resuhs in 



[ Vp,! yp,3 ] = H 



pA 



[ »7p,l ^P,3 ] 
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Taking the structure of this preamble into account and recalhng the assump- 
tion about the interference being mostly contributed from the first-order FT 

,-1 



neighbors, one can easily see that Cp i = ct, 

[ Vp,! Vp.S ] = -H'p,! 



[ ] 



HpiCpA2 



[ %,1 %.3 ] 



(18) 



where A2 is the orthogonal matrix 



and the pseudo-pilot Cp can be again pre-calculated as in the SISO case. An 
estimate of the CFR matrix at the subcarrier p can then be computed as 



H 



1 

2^ 



[ Vp,l Vp.,3 ] ^2 



(19) 



Remarks. 



1. It should be noted that, in all lAM variants except for lAM-I, Cp's have 
the same magnitude at all frequencies p. 

2. Note that, in view of the orthogonality of the matrix A2, noise enhance- 
ment is again controlled by the magnitude of the pseudo pilot Cp. Indeed, 

one can readily see that the covariance matrix of the noise term in ([H 
2 

equals j^-p/jy,; which corresponds to an MSB given by 



E \\H 



p,i 



■Hp J' 



N,ay\c„ 



(20) 



3. The preamble given above is easily generalized to systems with more than 
two transmit antennas, with A't being a power of 2. Simply select the 
corresponding matrix to be a Hadamard matrix of order iVt and 
build the preamble accordingly. The MSB will be still given by (|20p . 

4. For an A't x N,- system, the lAM preambles with guard symbols require 
2A't -I- 1 OFDM/OQAM symbols. This is equivalent to At + 1/2 com- 
plex OFDM/QAM symbols, that is, half a complex symbol longer than 
the preamble one would use in MIMO-OFDM. Nevertheless, even this ex- 
tra OQAM symbol can be avoided in practice, by eliminating the leading 
null symbol in the preamble above. This is often possible when interfer- 
ence from preceding data can be practically ignored, as for example in 
WiMAX [7], where the gap between two successive frames can play the 
role of the missing guard symbol. The number of complex OFDM sym- 
bols spent for the B-IAM-C preamble is somewhat larger: A't/2 more than 
that one would use for MIMO-OFDM. However, as it is demonstrated in 
the simulation examples, this increase in the preamble duration can offer 
considerable improvement in the estimation accuracy. 
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5. Eq. relies on the assumption that the prototype filter is sufhciently 
well localized in time so that interference from time instants other than 
the previous and next is negligible. In fact, often this is not the case in 
practice and there is nonnegligible interference to the {p, q) FT point from 
time instants g ± 2 as well. Hence the "pseudo pilot" matrix in (fT8|) is 
equal to 



instead, with corresponding to the ith transmit antenna and \Cp\ 7^ \Cp\. 
The latter implies that the two columns of the CFR matrix (i.e., the 
channels from the two transmit antennas) are not estimated with the 
same accuracy. The associated MSB for the subcarrier p is then equal to 

'^'^ (jjrp + J(?y) ' '^^^'^S^ ^'^^ analogous, yet somewhat more compli- 
cated, for Nt > 2. In general, the pseudo pilot matrix is no longer unitary. 
It only has (nearly) orthogonal rows. This fact was taken into account in 
the simulations; that is, the correct (not the ideal) pseudo pilot matrix was 
used in all cases. Note that this difficulty can be easily overcome, albeit 
at an extra cost in training bandwidth, if more than one guard symbol are 
placed among the A^t repetitions of the lAM preamble above. 



4.3 Optimal sparse preamble 

The MIMO preambles described above are of duration proportional to the num- 
ber of transmit antennas, iVt,. Shorter preambles, consisting of 3 OFDM/OQAM 
symbols per antenna (independently of iVt) as in the SISO case (with zero guard 
symbols surrounding a nonzero one), were also reported in [24 and are of the 
sparse type, i.e. with nulls among the pilot subcarriers. One of those pream- 
bles was based on earlier work on orthogonal training design for MIMO-OFDM 
systems [87] and shown to be also optimal in the MSB sense. The idea is briefly 
recalled here. Such a channel estimation approach was also recently adopted 
in [in]. 

Let N ^ he an integer]^ greater than or equal to 2A't0 Consider Nt 
sets of Lh pilots each, placed at subcarriers {pi,Pi + N,pi + 2N, . . . ,pi + {Lh — 
l)iV}, i — 1,2, . . . , Nt, with user-chosen starting positions pi,p2, ■ ■ ■ ,PNt G 
{0, 1, . . . , ~ 1} and with the rest of the subcarriers being nulled. Thus, the 
pilot subcarriers within each set are equispaced. Moreover, the pilot symbols 
are chosen to be equipowered. Pi's can be so selected as to have at least one 
zero between two nonzero pilots (sparse preamble), thus avoiding interference 
among them. 

Since the focus will be on the middle nonzero preamble vector in each an- 
tenna, the time index (normally q — 1) will be dropped for simplicity. Let 

n T 

T T T 
Vp Vp+N ' ' ■ yp+(L^_i)Ar , 

^If not, one can use for L/^ the smallest power of 2 that is not smaller than L^, that is, 
2Liog2 Lh\ . 

^^N > A'^t anyway, since the number of unknown parameters (NtNrLf^) cannot exceed the 
number of received signals (MNr). 
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with y,p = Up ^ defined as in (ITUl) , contain the received signals at the pth pilot- 
tone set and similarly define ry^^^ for the corresponding fi'equency domain noise 
components. Then one can readily verify that 



,(p) - 



ip) 



(21 



where 
with 



ip) 



(p) 



dia; 



iag (d'p 



with F{1 : N : M, 1 : Lh) being (in Matlab notation) the submatrix of the 
Afth-ordcr DFT matrix F consisting of its first Lh columns and every iVth of 
its rows. 



W'^P'^ =diag(l,< 



2p 



with wm — e^^~ denoting the i\/th root of unity, and h contains the im- 
pulse responses of the channels from transmit antenna i to all receive anten- 
nas. Then one can compute an estimate of the entire system impulse response 

h = [ ih-'Y {h-'Y 

system of equations 



• (/i-'^')^ ] by solving the NtN,Lh x NtN.Lh 



y(Pl) 




y(P2) 









*"PiVt 



Pi 

i2 



PNt 



Pi 



^PN, 





fj(pi) 


h + 


rq{P2) 






rj(p«t) 



y 



V 



or 



y^Ch + rj (22) 

Note that rj is white. Thus, the LS solution of ([^^ will be optimal in the MSE 
sense if C is unitary. A sufficient condition ensuring the latter is that 



di = d 



P+N 



= • • • = d; 



'p+{Lh-l)N 



, for all p, i 



(23) 



and the matrix 



D 



4. 



1^ 



4.. 



"PiVt 



be unitary. For example, with M — 8, Nt — 2 and Lf, 
Choosing D to be the orthogonal matrix 



(24) 

2, we have = 4. 



D 
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and pi = 0,p2 — 2, the middle (nonzero) vectors of the preambles employed at 
the two transmit antennas will be 



1 


1 


1 



1 





1 



-1 



1 



-1 





Note that, in the SISO case, the above reduces to the optimal sparse preamble- 
based approach of [49] as outlined previously. In fact, it can be verified that, in 
that case, the pilot symbols need only be equipowered (not as in (1231) 1 as proved 
in gg. 

5 Simulation Results 

In this section, simulation results are reported that demonstrate and compare 
the estimation performances of the lAM schemes previously discussed. In the 
SISO scenario, the iterative method of Hu et al. [IT] was also tested, with the 
preamble built as in Fig. (H^d), preceded by zeros (due, e.g., to the inter-frame 
gap) and followed by random i.i.d. data. Results of the POP method were not 
included since, as expected, it had a quite poor performance compared with the 
rest of the methods studied. For the sake of completeness, CP-OFDM is also 
included, with the minimum possible CP length (equal to the channel order). A 
realistic scenario where all the preambles are followed by pseudo-random dateF^ 
was considered. This is of importance for the results to be presented, since 
the assumption used above that the middle preamble OFDM/OQAM symbol 
receives negligible interference from the following data is not exact, even with 
the well localized prototype filters employed in the experiments. Thus, the late 
part of the preamble signal contains a portion of the front tail of the data burst 
as well, and this was taken into account when estimating the required transmit 
power. Note that this problem is not present in CP-OFDM and is due to the 
long tails of the OFDM/OQAM prototype filter impulse response. It should 
also be pointed out here that the SFB output signal due to the (SISO) 3-symbol 
OFDM/OQAM lAM preamble is of length {K + 1)M, which is approximately 
K + 1 times as long as that of CP-OFDM. One could consider shortening the 
preamble burst by eliminating its tails, as it was recently described in [T3]. As 
for the time invariance assumption for the channel and its validity throughout 
the preamble signal, one can see that, for scenarios similar to that of WiMAX 
for example [7], and for reasonable mobile speeds, the channel coherence interval 
spans many more than K + 1 M-blocks. 

The experiments utilized prototype filters designed as in [TT]. Filter banks 
with M = 512 and K = 3 were used. QPSK modulation was adopted. The 

^^In fact, what comes immediately after the preamble is usually control information rather 
than random data (e.g., [7]). Nevertheless, whatever structure is found in this control part of 
the frame is ignored here for the purposes of the experiments. 
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normalized MSE (NMSE), E H H where H is the CFR and 

H its estimate, is plotted with respect to the signal to noise ratio (SNR). The 
latter is defined as the ratio of the channel input (i.e., SFB output) power to 
the power of the noise. Note that in order to be fair with respect to the power 
transmission requirements, all the preambles are appropriately scaled so as to 
result in the same power at the SFB output. It must be emphasized that these 
power values turn out to differ much more with each other than they do at 
the SFB input. (In fact, earlier comparison studies seem not to have taken 
this into account and only equalized powers at the SFB input.) An example is 
shown in Fig. [HI The figure also shows that, as pointed out in [SD], this kind 
of lAM preambles result in high PAPR signals at the SFB output, due to their 
deterministic/periodic nature. In fact, as it can be seen in the example, this 
effect is more evident in lAM-C than in lAM-R, while it is much less severe in 
lAM-I due to its pseudorandom symbols. E-IAM-C exhibits even higher peaks 
than lAM-C. These peaks were included in the average power estimates of the 
corresponding channel inputs. 



5.1 SISO Systems 

Training in CP-OFDM was based on a single-symbol complex QPSK pseudo- 
random preamble. The iterations in the method of |41j relied on a single data 
OFDM/OQAM symbol following the pilot symbol. At moderate to high SNR 
values, two iterations sufficed for convergence. Figs.lTJa) and (b) show the MSE 
performance of the methods under study for Veh-A and Veh-B channels |75| , re- 
spectively. Observe that in both cases E-IAM-C outperforms the other schemes 
in the whole SNR range considered, with the iterative method of Hu et al. ex- 
hibiting the worst performance among the OFDM/OQAM schemes. Moreover, 
all the OFDM/OQAM methods perform better than CP-OFDM at low to mod- 
erate SNR values. At higher SNRs, CP-OFDM takes over while the NMSE 
curves of the OFDM/OQAM schemes reach an error floor. This is a well known 
phenomenon in such systems and is due to the fact that the approximation ([5]) 
is not exact for channels of significant time dispersion. Hence there is residual 
intrinsic interference, which is hidden by the noise at low SNRs and shows up 
in the weak noise regime. It should also be noticed that the OFDM/OQAM 
methods are less well performing with channels of higher frequency selectivity 
and the crossing point with the CP-OFDM curve appears earlier in that case 
(see Fig. [2llb)). This behavior is again explained by the fact that, in the highly 
frequency selective case, the assumption of a locally flat CFR, and hence ([5]), is 
much less valid. 

Fig. [5] provides an example of the performances obtained when the pream- 
bles are only normalized to equal powers at the SFB input (instead of at the 
transmitted signal as above) , as it is the case in most of the existing literature. 
Comparing with Fig. [TJ^a), one can see that the performance of all methods is 
then improved (with the error floor effect being significantly remedied) while 
the differences among them are reduced. 
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5.2 MIMO Systems 



The channel taps were modeled in accordance with the Kronecker model, that 
is h{l) = RI^'^ {l)h^(l) (l) , I = 0,1, . . . , Lfi — I, where the receive and trans- 
mit correlation matrices Rt{l), R^i^) ^-re generated according to the exponential 
model, and h^{l) is a matrix of i.i.d. zero mean circularly symmetric Gaussian 
entries following a Veh-A or Veh-B power delay profile. Weak spatial correla- 
tion was assumed at both the transmit and receive sides of the channel. The 
CP-OFDM preamble was similarly chosen, namely Ajy (X) x with x being a 
pseudorandom Af -vector of complex QPSK pilots. 

Fig. [S] shows the MSE performance of the methods under study in 2 x 2 
systems. Veh-A and Veh-B channels were considered in Figs. [SJa) and (b), 
respectively. Once more, E-IAM-C outperforms the other schemes in the whole 
SNR range considered. Moreover, all the lAM methods perform better than CP- 
OFDM at low to moderate SNR values. At higher SNRs, the error floor effect, 
which was observed in the SISO scenario, appears here too. As previously, it is 
more severe for more frequency selective channels (cf. Fig. IHl^b)). 

6 Conclusions 

The problem of preamble-based channel estimation in OFDM/OQAM systems 
was revisited in this paper, and a review of the preamble structures and associ- 
ated estimation methods was given, for both SISO and MIMO systems. Results 
on MSE-optimal design of preambles with zero guard symbols were also re- 
viewed for SISO systems. Analogous results for the MIMO scenario are still to 
be derived. The estimation performances of the principal methods (including 
lAM variants and the iterative method of [H] ) were demonstrated via results of 
simulation in a realistic (including interference from data) and fair scenario and 
for both mildly and highly frequency-selective channels. A recently developed 
method of the lAM type, E-IAM-C, was seen to outperform the rest under such 
conditions and at all considered SNR values. Further investigation is required 
to assess the loss in performance from tail truncation in the lAM preambles 
and mitigate the high PAPR problem associated with their periodic structure. 
Moreover, improved I AM designs as well as optimality results are needed for 
MIMO-OFDM/OQAM. 
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A Time- Frequency Neighborhood 

In this appendix, the interference terms from the neighbors of a FT point are 
computed in closed form, based on the knowledge of the prototype filter. More- 
over, their symmetries are demonstrated. 

These are the quantities we have previously denoted as 

{9} m-\-p,n+q ^{{9ni+p,7i+q\g7n,n)} 

where 

{9m+p,n+q\gm,n) — Sm+p.n+g (Offm,n (0 ■ 
I 

The latter has to be computed for p G {—2, —1, 0, 1, 2} and q <E {—1, 0, 1}, with 
{p, q) 7^ (0, 0). In general, it can be written as: 

{gm+p,n+q\9m,n) = X! 5' " (^^ + 9)"^^ 9 ~ e^~^(' ^) e^'^'^ EE J- (5) "f^ 

(25) 

where Aip = ipm+p,n+q — ^m,n = (p + " {mq + pn + pq)Tr . Substituting this 
in psp and after some algebra, the following simplified form results: 

Lg-l-max(gf-.O) 

{g„.+p,n+q\9m,n) = (-1)"^'^ e'^^^^ ^1 9(1)9 [1 + 1y) 

!=inax(-g^,0) 

(26) 

Observe that this is independent of n. 

Let us consider the various cases in detail: 
l)q = 0,p€{±l,±2} 

With p — ±1, (pS)) becomes 

1=0 

and hence 

(5)™£,„ = ±6^^^^ 9\l)e^'^^'- (27) 

1=0 

In view of the fact that the latter quantities are real and hence equal to their 
conjugates, we can write 

/3^(ff)™+"i,„ = -(3)™L"i,„- (28) 

For p — ±2, we have 

{gm±2,n\gm,n) = -e^^^f^-'-i' ^ g'{l)e^^^'' = -eT^^(^-'-i)G(2)(±2), 

1=0 

with G^^''{p) denoting the Af -point discrete Fourier transform of g. This coin- 
cides with the ambiguity function of 5, Ag{q,p), calculated at (0,±2). But the 



20 



latter is known to equal zero for all g's that can be used as prototype filters in 
such filter banks [3S]. Hence, 

/ \m.n n 
(5')m±2,n = 

2) <Z = -l,pe {-2,-1,0,1,2} 
With (J = -1, (IMl) becomes 

(ff™+p,„-i|5™,«) = E -9(0.9 - f ) e^^^'. 



Thus, 

io-l 



l=M/2 ^ 

= -(-ir7- (29) 



For p = ±1: 



2 , 

l=M/2 ^ ' 

(-1)"^, (30) 



while for p = ±2: 



i=Af/2 ^ ^ 

= (-l)"e. (31) 
3) g= l,pe {-2,-1,0,1,2} 
Substituting q = 1 in ([^ yields 

(ffm+p.«+i|ffm,„) = J • /(-l)'"e-^f^^ E ' 5(05 + y) eJ^-^' 
and after a change of variable in the summation: 



(ff™+p.„+i|3™,„) - J • /(-l)™+J'e-^f^^ E 5(05 - y) 



Thus, 



2 



(.9)"±"i,„+i = -,7±^(-l)'"eT^^"^ E 5(03(^-^)e±^^' 



'£5(05(^-1) 



l=M/2 

La-1 



(5)™f2,„+i = -(-l)-eT^^(--^' E .9(05 (/-y)^^''^' 

i=Af/2 
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which, in view of the above, result in: 



(-I)'"!, (32) 
i-irS, (33) 
-(-l)'"e. (34) 

It thus turns out that the 5x3 neighborhood of the (m, n) point is as follows: 



i-ire 





-(- 


-l)'"e 


(-1)'"(5 


-/3 


(- 


-1)"'6 


-(-1)^7 




(- 


■1)'"7 


(-1)™^ 
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(- 


-1)"'S 


(-l)™e 





-(- 


-l)"^e 
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Figure 2: The MIMO-OFDM/OQAM system. 
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Figure 3: Preamble structures for (a) lAM-R , (b) lAM-I, (c) lAM-C methods, 
and E-IAM-C method for (d) e > and (e) e < 0. M = 8. OQPSK modulation 
is assumed, do, di are randomly chosen BPSK symbols. 
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Figure 4: Preamble structures aiming at interference cancellation: (a) [H]; (b) 
SS]; (c) [5S], and avoidance: (d) gH US]. M = 8. OQPSK modulation is 
assumed. 
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Figure 5: lAM-C preamble for a 2 x X system, with (a) and (b) correspondig 
to the two transmit antennas. M = 8. OQPSK modulation is assumed. 
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Figure 6: Magnitude squared of the modulated lAM preambles. AI — 512, 
K = 3. QPSK modulation. 
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Figure 7: Estimation (NMSE) performance of the preamble-based methods for 
SISO systems. Filter banks with M = 512, K — 3 have been used. Rayleigh 
channels were assumed, with a (a) Veh-A and (b) Veh-B profile. 
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Figure 8: As in Fig. [Tl^a), but with preamble powers equalized only at the SFB 
input. 
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Figure 9: Estimation (NMSE) performance of the I AM preamble-based methods 
for MIMO 2x2 systems. CP-OFDM is included for the sake of comparison. 
Filter banks with M = 512, A' = 3 have been used. Rayleigh channels of (a) 
Veh-A and (b) Veh-B profiles have been assumed. 
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